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ABSTRACT 

We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts mea- 
sured from Sloan Digital Sky Survey III Data Release Eight (SDSS DR8). The galaxies have photo- 
metric redshifts between z = 0.45 and z — 0.65, and cover 10,000 square degrees and thus probe a 
volume of 3/i~ 3 Gpc 3 , enabling tight constraints to be derived on the amount of dark matter in the 
form of massive neutrinos. A new bound on the sum of neutrino masses < 0.26 eV, at 95% 

confidence level (CL), is obtained after combining our sample of galaxies, which we call "CMASS", 
with WMAP 7 year Cosmic Microwave Background (CMB) data and the most recent measurement 
of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a 
conservative multipole range choice of 30 < I < 200 in order to minimize non-linearities, and a free 
bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy 
bias model using mock catalogs, and find that this model causes a small (~ 1 — 1.5c) bias in f^DM^ 2 - 
For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing 
additional, shot noise-like free parameters. In this conservative case, the bounds are significantly 
weakened, e.g. ^m„< 0.36 eV (95% confidence level) for WMAP+HST+CMASS (£ max = 200). We 
also study the dependence of the neutrino bound on multipole range (£ max = 150 vs £ meiK = 200) and 
on which combination of data sets is included as a prior. The addition of supernova and/or Baryon 
Acoustic Oscillation data does not significantly improve the neutrino mass bound once the HST prior 
is included. 



A companion paper (Ho et al. 2012) describes the construction of the angular power spectra in 
detail and derives constraints on a general cosmological model, including the dark energy equation of 



state w and the spatial curvature fi#, while a second companion paper Seo et al. ( |2012 ) presents a 
measurement of the scal e of baryon acous tic oscillations from the same data set. All three works are 
based on the catalog by Ross et al. (2011). 



1. INTRODUCTION 

During the last several years, experiments involving so- 
lar, atmospheric, reactor and accelerator neutrinos have 
adduced robust evidence for flavor change, implying non- 
zero n eutrino mass, see Ref. Gonzalez-Garcia fc Maltoni 
( 2008 ) and references therein. The most economical de- 
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scription of the neutrino oscillation phenomena requires 
at least two massive neutrino mass eigenstates to explain 
the two mass differenc e!^ Am?, = 7 . 59 ■ 1 ~ 5 eV 2 and 
Am|a = 2.5 ■ 10~ 3 eV 2 |Schwetz et al.| ( |201l| ; |Fogli et al . 



(2011), which drive the solar and atmospheric transi 



tions. Despite the remarkable success of past and present 
oscillation experiments, and the promising prospects for 
future searches, the individual neutrino masses and the 
Dirac versus Majorana neutrino character are key ques- 
tions that continue to be unanswered by oscillation ex- 
periments. 

Direct information on the absolute scale of neutrino 
masses can be extracted from kinematical studies of weak 
decays producing neutrinos. The present upper bound on 
the electron-neutrino mass from tritium beta -decay ex- 
perim e nts is 2 eV ( 95% confidence level (CL)) Lobashev 
( |2003| ; |Eitei| ( |2005j ) , and in the future the KAl'RIN ex- 
periment is expected to be sensitive to e lectron-neutrino 
masses appro aching 0.2 eV (90% CL) Otten & Wein- 
heimer (2008). Searches for the Majorana neutrino na- 
ture involve neutrinoless double beta decay f3f3(Qv), a 
rare and as yet unobserved transition between two nu- 

13 Neutrino oscillations are described by mass squared differ- 
ences and not by the absolute values of the mass eigenstates 



2 



clei. Observational upper limits on (3[3(0v) rates pro- 
vide an upper bound on the so-called effective Majorana 
mass of the electron neutrino, (m e s) < 0.3 — 1.0 eV, 
bound which would only apply if neutrino s are Majo- 
rana particles Gomez-Cadenas et al. (2011). Forthcom- 



ing (3f3(0v) experiments aim for sensitivity approaching 
(m c ff) < 0.05 eV Gomez-Cadenas et al. (2011). 

Cosmology provides one of the means to tackle the ab- 
solute scale of neutrino masses. Some of the earliest cos- 
mological bounds on neutrino masses followed from the 
requirement that massive relic neutrinos, present today 
in the expected numbers, do not saturate the critical den- 
sity of the Universe, i.e., that the neutrino mass energy 
density given by 



93.1/i 2 eV 



(1) 



satisfies Vl u < 1. The Universe therefore offers a new lab- 
oratory for testing neutrino masses and neutrino physics. 
Accurate measurements of the Cosmic Microwave Back- 
ground (CMB) temperature and polarization anisotropy 
from satellite, balloon-borne and ground-based experi- 
ments have fully confirmed the predictions of the stan- 
dar d cosmological model and all ow us to weigh neutri- 



nos [Lesgomgues&Pastor (2006). Indeed, neutrinos can 
play a relevant role in large-scale structure formation 
and leave key signatures in several cosmological data 
sets. More specifically, the amount of primordial rela- 
tivistic neutrinos changes the epoch of matter-radiation 
equality, leaving an imprint on CMB anisotropics. Af- 
ter becoming non-relativistic, their free-streaming nature 
damps power on small scales, suppressing the growth 
of matter density fluctuations and thus affecting both 
the CMB and g alaxy clustering observables in the low- 
redshift universe |Lesgourgues fc Pasto r (2006). Measure- 
ments of all of these observations have been used to place 
new b o unds o n neutrino phys i cs from cosmology Elgaroy 
eTal] (|2002 ]) ; |Spergel etaL] ([2003) ) ; I Hannestadl (|2003J; 
Allen et al.| ([2003 ): Tegmark eTaT l (|2004D ; IBarger et al. 



(2004); Hannestad & Raifelt (2004); Crotty et al. (2004) 
Seljak et all (|2005|); [Elgaroy & Lahavl j2005|r | Hannes 
tad| (|2005|);]G ooDaj et ah] J2006I); |S^geTetal '| (| 2007[) 
Seljak et al.| (|2006[) ; |Fogfi est al.[(|2008p ; IKomatsu et al 



'([2009 [) ; |Rei(j eTal ] (|M0 ] ); JReid' e t ^'|2V)10a[) 



(M~iiTT(|2TTTI7C:|Keid efaITi|2010bD; IKo matsu et al 
Sait o" eT al.| (|201 1[); |Kiemer-S0rensen et al.| (12011 



1' homa sl 
|(|2011|); 



Bcn- 



son et al.| ( |201ip , with a current limit on the sum of neu- 
trino m asses Y>m u ^ 0.6 eV at 95% CL (e.g. Reid et al. 



(2010)), depending on the precise combination of data 



sets and on the underlying cosmological model. 

We present here neutrino mass bounds from the fi- 
nal imagin g data set o f the S loan Digital Sky Survey 



(SDSS-III) York et al. (2000) using the photo metric 



redshift catalog of Hoss et al. Ross et al.| J20 11). We 
consider the CMASS sample Whit e~eTlilT(|2TIllll o f lu- 
minous galaxies of SDSS DRij |Aihara et al.| ( |20ll[ ), the 
eighth data release of SDSS and the first data release 
of the Baryon O s cillati on Spectroscopic Survey (BOSS) 



(20111, with photometric redshifts 



Eisenstein et al 
0.45 — 0.65. This sample covers an area of approximately 
10,000 square degrees and consists of 900,000 galaxies. 
It is thus the largest sample of luminous galaxies so far 
and promises strong constraints on the neutrino proper- 



ties (see Thomas et al. (2010) for an analysis of a slightly 
smaller SDSS photometric sample). 

We derive neutrino constraints from the angular power 
spectra of the galaxy density at different redshifts, in 
combination with priors from the CMB and from mea- 
surements of the Hubble parameter, supernovae dis- 
tances and the BAO scale. The spectra and the analysis 
of a minimal ACDM c osmology are de scribed in detail in 



our companion paper 
ment of the BAO sea 



Ho et al. (2012) and the measure- 
e from the spectra is pr esented in 



Seo et al. (2012). We will 



a separate companion paper 

often refer to these works for details and focus here on 
the neutrino bound. 

The structure of the paper is as follows. In section [2[ 
we describe the data set and the derived angular spectra. 
We then discuss our theoretical model for the spectra and 
their cosmology dependence in section [3] In section [4] we 
explain the specific signature of neutrino mass on galaxy 
clustering data. We test our model for the angular power 
spectra against mock data in section [5] and present the 
constraints on the sum of the neutrino masses and other 
parameters for several data combinations in section [6[ 
Finally, we discuss these results and conclude in section 

m 

2. DATA 

The data and the method for obtain ing angular spectra 
have b een descr i bed in detail in Ref. Ross et al. (20111 
and in Ho et al. (2012 1. Here we provide a brief descrip- 
tion of the main properties and refer the reader to those 
papers for details. Our galaxy sample i s obta ined from 
im aging data from D R8 Aihara et al[ ( |2011[ ) of SDSS- 
III j^keFal]([2000]). This survey mappecTabout 15,000 
square degrees of t he sky in five pass bands (u, g, r, i and 
z) iFukugita et al. ( 1996| ) using a wide field CCD camera 
Uunn et al. ( |1998 ) on the 2.5 meter Sloan te lescope at 



Apache Point Observatory Gunn et al. ( 2006| (the sub- 
sequent astr qmetric calibration of these imaging data is 
described in |Pier et al.J (|2003 [)). A sample of 11 2, 778 
galaxy spectra from BOSS" Eisenstein et al. (20111 were 
used as a training sam ple for the photo metric redshift 
catalog, as described in Ross et al. (2011). 
We focus on the approximately stellar mass- limited 



CMASS sam ple of luminous galaxies, detailed in White 



et al. (2011), which are divided into four photometric 
recTsMfih5ms, z photo = 0.45 - 0.5 - 0.55 - 0.6 - 0.65. 
The photometric redshift error lies in the range a z (z) — 
0.04 — 0.06, increasing from low to high redshift. Figure 
Q] shows the estimated true redshift distribution of each 
bin, d etermined using t he methods described in section 



5.3 of Ross et al. (2011) 



The calculation of the angular power spectru m for each 
bin is described in detail in our companion paper |Ho et al.| 
(20121 and uses th e optimal qua d ratic estimatorT OQE) 



meThod outlined in Scljak (1998); Tegmark et al. (1998) 



Padmanabhan et alT 1 2003[ 2007). The four power spec 
tra are binned in I space with a typical wave band width 
of A£ = 10. The expectation value of the power spec- 
trum in a wave band is a convolution of the true power 
spectrum with a window function of width roughly equal 
to the typical wave band width. Examples of these win- 
dow functions a re shown in Fig. 3 of the companion paper 

r 



Seo et al 
theoretical 



[2012 ). When fitting the data to the underlying 
model, we always apply these window func- 
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Fig. 1. — Normalized true rcdshift distribution of CMASS galaxies in four photometric redshift bins. The number of galaxies in each bin 
is 214971, 258736, 248895 and 150319 (from low to high redshift). 



tions to the theoretical power spectra before calculating 
the likelihood relative to the data. 

The four power spectra are plotted with their error bars 
in Fig. [2j The solid curves represent theoretical power 
spectra based on several assumptions. These power spec- 
tra will be discussed in detail in section [3j Sinc e the 
low I wave bands are more prone to systematics |Ross 



etal 
bands 



(2011), we are conservative and do not consider 
with 



< 30 in our analysis. We shall apply cuts 
at ^ max = 150 and 200 in order to suppress uncertain- 
ties from non-linear corrections to the modeled power 
spectra, as discussed in the following sections. The me- 
dian redshift (z ~ 0.55) contributions to these maximum 
angular modes arise from three-dimensional modes with 
wave vectors k « O.lO/iMpc^ 1 and k « 0.14/iMpc -1 , re- 



spectively, 
slice for £ ma . 



We thus use 17 (12) 
= 200(150). 



data points per redshift 



3. MODELING THE ANGULAR POWER SPECTRA 

The galaxy overdensity in the i-th redshift bin can be 
expanded in terms of spherical harmonics, 



<5«(n) 



t in 



(2) 



so that the angular power (and cross) spectra are defined 
as 



,(•) „(i> 



C, 



(«) riK cK 

°U> °mm' ' 



(3) 



where 6$ is the Kronecker delta function. As mentioned 
in the previous section, we do not estimate our spectra by 
directly transforming the observed density field to har- 
monics space, but use instead the optimal quadratic es- 
timator technique. To constrain the sum of the neutrino 
masses and other cosmological parameters, the observed 
spectra are compared to a cosmology dependent model, 
which we now describe. 

The galaxy overdensity on the sky is a line-of-sight 
projection of the three-dimensional redshift space galaxy 
overdensity 8 g {d(z) n, z), 



where 



«W(ft)= I dz 9i (z) (6 g (d(z)h 1 z)-(H(z))- 1 h-V(n-v(d(z)h,z)) 

(4) 

drii/dz(z) 



J dz' drii/dz(z') 



(5) 



is the normalized redshift distribution of galaxies in bin i 
(with drii/dz{z) the number of galaxies per steradian per 



4 




HF 

linear 

Limber (no RSD) 
no RSD 
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linear 

Limber (no RSD) 
no RSD 




6 4 =2.40, WMAP7+f/„ best-fit cosmology 
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Fig. 2. — Observed power spectra (black points) with error bars and theoretical power spectra (solid curves). We show the theoretical 
power spectra for different models: the default, HaloFit (HF) based model used in our analysis (black; see text for details), the same 
model, but using the linear matter power spectrum as input (red), the default model, but using the Limber approximation (blue) and the 
default model without rcdshift space distortions (green). We restrict ourselves to the range i = 30 — 200 in our analysis. For the theoretical 
spectra, we assume the WMAP7+HST best-fit cosmology and use the bias bj that best fits the data. We do not here include the shot noise 
parameters a^. 



unit redshift), d(z) is the comoving distance to redshift 
z (assuming a flat universe) and v is the galaxy velocity 
field. The velocity term arises because gradients of the 
peculiar velocity contribution to the distance in redshift 
space change the volume, and consequently, the number 
densit3p^| 

We assume a linear, scale-independent bias for the 
galaxy density, 



where f3(z) = f(z)/b g (z) is the redshift distortion pa- 
rameter and 



_ dlnD(z) 
}[Z)= d\na 



(8) 



5 g (x,z) = b g (z)S m (x,z) 



(6) 



is the growth factor (with D(z) the linear growth func- 
tion). In the presence of neutrinos, the growth function 
is no longer scale independent at late time as the neutri- 
nos su ppress growth on scales below the free streamin g 
length |Hu & Eisenstein| (|1998|) ; |Eisenstein & Hu| ([1997} 



with 8 m the matter overdensity. For the peculiar velocity 
field, we use the continuity equation in the linear regime, 
which gives for a Fourier mode with wave vector k, 



-i/3(z)6 g (k) ¥ , 



(7) 



14 Instead of writing the projected galaxy overdensity as an inte- 
gral over the observed redshift (including peculiar velocity contri- 
butions) as in Eq. 14b, one could equivale ntly do the in t egral over 
"true" cosmic redsmft, see Ref. |Padmanabhan et al.| | |2007[ ), in 
which case only the true three-dimensional galaxy overdensity ap- 
pears explicitly and the redshift space distortions come in through 
a modification of the distribution gi(z). 



but not on larger scales (with a broad transition regime 
in between) . We shall ignore the scale dependent growth 
in j3(z) since it is a small (<?C 10%) correction to the 
already small effect (on the scales of interest here) of 
redshift space distortions (RSD, see Fig. [2]). However, 
this scale-dependent growth is included in the real space 
power spectrum, as this is the main signature of massive 
neutrinos. 

We simplify our treatment of the galaxy bias by follow- 
ing the approach of our two companion papers, adding 
four free parameters bi to describe the bias in each bin. 
The results from our simulations barely change when 
considering a bias b g {z) that varies within redshift bins, 



■5 



showing that this is a sa fe approxima t iorpl It then fol- 



lows from the above (see Fisher et al. | rtl994| i ; jHeavens fc 
Taylor| ( |1995[ ); |Padmanabhan et al.| p007| j thaT~ 



C 



b 2 ~[ k 2 dkP m (k,z = 0) (Af\k) + Af SD ^(k)) 2 , 

(9) 

where P m (k,z = 0) is the matter power spectrum at 
redshift zero and 

Af\k) = J dz 9i (z)T(k,z)M(kd(z)) . (10) 

Here, ji is the spherical Bessel function and T(k, z) the 
matter transfer function relative to redshift zer cED The 
contribution due to redshift space distortions is 



A RSD,(i) (fe)=ft / dzgi{z)T{kjZ) 



1(1-1) 

(2Z-1)(2Z+1) 

(i + m + 2) 

(2l + l)(2l + 3) 



{2l 2 + 21 - 1) 
(2/ + 3) (2/- 1) 

Jl-2(kd(z)) 



ji(kd(z)) 



ji +2 (kd(z)) 



(11 



To compute the matter power spectrum at a given red- 
shift P m (k, z) = P m (k, z = 0) T 2 (k, z), we first make use 



of the CAMB code Lewis et al. (2000), which provides 
the linear power spectrum by integrating the Boltzmann 
equations of all species including massive neutrinos. We 



then apply the HaloFit prescription^] Smith et al. ( 2003 1 
to the linear power spectrum to account tor non-linear et- 
fects on the matter power spectrum. 

While in the linear regime the galaxy spectrum is easy 
to model, calculations on non-linear scales inevitably 
have large uncertainties. This effect is aggravated by the 
presence of massive neutrinos since for the massive neu- 
trino case the non-linear regime has been explored less 
extensively in the literature than for a vanilla ACDM 
model. In the non-linear regime, the matter power spec- 
trum receives corrections due to gravitational collapse, 
the galaxy bias becomes scale-dependent, and redshift 
space distortions receive important contributions from 
velocity dispersion in collapsed objects. We take into 
account non-linear corrections to the matter spectrum 
using HaloFit. The effect of non-linearities on redshift 
space distortions at the relevant scales here is small as 
it is largely washed out by line-of-sight projection. How- 
ever, we do expect significant corrections to our model 
on small scales due to non-linear galaxy bias, which we 
address below. 

For angular scales where non-linear effects cannot be 
ignored, the contribution to a given angular mode £ from 
a redshift z comes exclusively from three-dimensional 
modes with wave vector k w £/d(z). Clearly, to avoid 

15 A similar approach is considered to model /3(z), appearing in 
the redshift space distortion contribution. For each bin we calculate 
an effective growth rate fi = (^DM^i)) ' 56 where Zi is the mean 
redshift of the i-th bin, ignoring the scale dependence of the RSD 
growth. 

16 The transfer function is defined as 5 m (k, z) = 
T(k,z) S m (k, z = 0). 

17 Recently, |Bird et aI.|j20lT| ) developed an extension to HaloFit 
that incorporates the eitect ot massive neutrinos. We do not use 
this prescription as the correction to standard HaloFit is negligible 
on the scales of our interest. 



large non-linear corrections, the analysis must be re- 
stricted to low £. On the other hand, the density of 
modes per unit £ increases with £ so we want to use as 
many modes as possible without biasing the results. Fig- 
ure [3] (left panel) depicts (as a function of redshift z) the 
value of £ above which non-linear corrections to the three- 
dimensional power spectrum contributions to the angu- 
lar spectrum become important (i.e. £nl = ^nl(^) d(z)), 
considering various assumptions for the non-linear scale 
^nl(-z)- Given that most of our signal is produced in the 
range z — 0.45 — 0.65, and assuming that our model be- 
comes inadequate at k > 0.15/iMpc -1 , we conclude that 
a conservative choice would be £ max somewhere between 
150 and 200. 

Alternatively, we can obtain an indication of the im- 
portance of non-linear galaxy bias by considering the ef- 
fect of non-linear corrections to the matter power spec- 
trunpj (which we do include in our model) . The right 
panelm Figure [3] therefore shows the signal to noise ra- 
tio squared in the difference between our default model 
and the same model, but using the linear matter power 
xSpectrum instead of the non-linear (HaloFit) one. The 
'signal to noise reaches one somewhere between £ max = 
150 and 200, corresponding to contributions of modes 
fcmax ~ 0.10/iMpc -1 and fc max sa 0.14/iMpc" 1 at the 
median redshift z — 0.55. Finally, a more concrete in- 
dication of the importance of non-linear galaxy bias to 
the rang e of scales of our cho ice is given by Fig. 13 (left 
panel) of |Hamaus et al.| ( |2010[) , which shows the halo bias 
as a function ot three-dimensional mode k. Since for our 
sample of galaxies the bias b ~ 2, the plot confirms that 
there is only a mild bias variation in the relevant range of 
three-dimensional scales relevant to the multipole range 
we have chosen. 

Based on the above discussion, we choose a default 
value £ max = 200, but we will also present results for 
the more conservative choice f ma x = 150. While it is 
possible to model the gal axy spectrum in a more sophis- 
ticated manner (see e.g. |Saito et al.] ( |2008[ |2009[ |2011 1 
for an approach based on perturb ation theory and the 
local b ias model iMcDonaldl ( 2006 1 , and Swanson et al. 



(20101 for a cross-comparison ot a number ot methods) 



we consider it appropriate, given the multipole range we 
include, to use the simple model described in Eq. Q, 
characterized by bias parameters hi. In addition to this 
model, we also consider an alternative model with more 
freedom, by adding shot noise-like parameters at (see also 
our companion papers), 



(»») 



k 2 dkP m (k,z = 0) (fc) + A 



RSD,(i) 



(k) 

(12) 

The parameters serve to mimic effects of scale- 
dependent galaxy bias and to model the effect of poten- 
tial insufficient shot noise subtraction. This model is a 
versi on of what is sometim e s referred to as the "P -model" 
(e.g. Hamann et al. (2008); Swanson et al. ( 2010 )) and is 
indep e ndently motivated by th e halo m odel Seljak (2000 
20011); iSchulz fc Wh itel fl2006|); |Guzik et a l.| (|2007T) anc 



the local bias ansatz |Scherrer fe Weinberg (1998); Coles| 



18 However, one must keep in mind that this may underestimate 
the eitect ol non-linear galaxy bias, as' galaxics arc more st rongly 
clustered than matter and arc thus prone to larger non-linear cor- 



rections. 
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Fig. 3. — Left panel: Minimum multipolc at which 3-D power spectrum contribution to the angular power spectrum receives important 
non-linear corrections, as a junction of redshift, Inl = ku^d^^^Nc consider several choices of the non-linear scale /cnl ■ The dashed curves 



0)/Rnl.( z ) x 0-lfeMp c 1 , where -Rnl( 2 ) is such that the m atter overdensity variance averag ed over spheres with this radius equals one 
(using the linear powe r spectrum). Right panel: The di fference as a functio n of l max between our default template, which u ses Halofit 
and a template using the linear matter power spectrum, given the covariancc matrix tor the (JMASS spectra. We assume the WMAF7 plus 
HST best fit cosmology and fix the bias parameters 6, = 2 (a; = 0). Both plots suggest that non-linear effects start to become (mildly) 
relevant at l m ax between tbU and 2UU. 



et al | ([1999b; |Saito et al.| <|2009}. We further discuss the 
validity of the parameterizations with and without a, in 
section [5] 

Figure [2] shows the theoretical galaxy spectra as de- 
scribed in this section. The error bars follow from the 
optimal quadratic est imator method used to construct 
the power spectra (see Ho et al. (2012 1 for details). Com- 
paring the spectra with (black) and without (green) red- 
shift space distortions shows that this effect is negligible 
for t > 50 and is probably not relevant for the range of 
scales we use in our data analysis, i.e. I > 30. Although 
we never employ it, we also show the effect of using the 
Limber approximation |Limber (1954) and find that for 
£ > 30 it works excellently. 

4. COSMOLOGICAL SIGNATURE OF NEUTRINOS 

In the analysis presented in this paper we assume that 
there are three species of massive neutrinos with equal 
masses m„. Massive neutrinos affect galaxy formation at 
scales below the Hubble horizon when they become non 
relativistic, 



can not grow due to the large neutrino velocity disper- 
sion. As a consequence, the growth rate of density per- 
turbations decreases and the matter power spectrum is 
suppressed at k > kf s . Since the free streaming scale de- 
pends on the individual neutrino masses and not on their 
sum, a measurement of kf s could, in principle, provide 
the ordering of the neutrino mass spectrum. In practice, 
such a task appears extremely challenging, see |Jimenez] 



0.0145 



f^DM WVlpC 



(13) 



with f^DM the present total dark matter energy density, 
i.e. cold dark matter plus massive neutrinos, relative to 
the critical density. The non-relativistic neutrino over- 
densities cluster at a given redshift z only at scales where 
the wavenumber of perturbations is below the neutrino 
free streaming scale 



k fs (z) 



0.677 



(l + z)V2 Vl eV 



\A^DM /lMpC~ 



(14) 



due to the pressure gradient, which prevents gravita- 
tional clustering. On spatial scales larger than the free 
streaming scale k < kf s , neutrinos cluster and behave as 
cold dark matter (and baryons). Perturbations with co- 
moving wavenumber larger than the free streaming scale 



et ah] ([2010 1. 

Figure E] illustrates the effect of massive neutrinos on 
the angular power spectra. The solid curves depict 
the results for the four redshift bins exploited here in 
the case of a ACDM model assuming no massive neu- 
trino sp_ed£s_and^est_iit parameters to WMAP7 year 
data Larson et al. (|2010 ' )• Komatsu et al. (2010) and 
HST-77B data |Riess et al.T ( |2vJll| ). The dashed curves de- 
note the angular power spectra results assuming three 
massive neutrinos with ^ m u = 0.3 eV and keeping the 
cold dark matter mass energy density constant. In the 
presence of massive neutrinos the angular power spectra 
are suppressed at each redshift at an angular scale that 
is related to the free streaming scale by t ~ d(z) kf s (z). 
Therefore, the larger the neutrino mass (or the redshift), 
the larger the lowest angular wavenumber at which the 
power spectrum is maximally suppressed. In the redshift 
range of interest here and for ^2 m v — 0.3 eV, the sup- 
pression angular scale appears in the range I = 20 — 50 
(however, there is some suppression even at lower £). 
Note as well that there will exist a strong degeneracy be- 
tween the neutrino masses and the amount of cold dark 
matter, since, in principle, one could partially compen- 
sate the growth suppression induced by massive neutri- 
nos at scales k > kf s by increasing the cold dark matter 
mass-energy density. Combination with CMB and H 
data will help to break this degeneracy. 

Neutrino masses affect the angular power spectra C^, 
see Eq. in two different ways: suppressing galaxy 
clustering and the growth of structure via P m (k) as well 
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as modifying the background expansion rate via the co- 
moving distance which appears in the argument of the 
Bessel function jg. Among these two effects (i.e. growth 
versus background) we find that the growth suppression 
effects in the matter power spectrum due to the presence 
of massive neutrinos will dominate over background ef- 
fects. Therefore, the neutrino mass constraints presented 
in the following analysis arise mostly from the suppres- 
sion of clustering rather than from purely geometrical 
effects. 

5. MOCKS 

We first consider angular spectra based on mock galaxy 
catalogs to test that neither our method of estimating 
the spectra nor our modeling of the spectra introduces a 
bias in the reconstructed cosmological parameters. We 
use twenty independent CMASS mock catalogs based on 
N-body simulations and a Halo Occupat i on Di stribution 
(HOD) model describ ed in White et al.j (|2"011[) (see also 



Hoet al.| ( |20f2| r r[Seo et al. 



our companion papers 

for details). The input cosmology for the Simula 
fl DM = 0.274, h = 0.7, n b = 0.046, n s = 0.95, a 8 



(2012) 



ions is 
= 0.8 

in a spatially flat universe, with n s and ag the scalar 
spectral index and the linear rms density fluctuations in 
spheres of radius 8ft. _1 Mpc at z = 0, respectively. Neu- 
trinos are massless in the input cosmology. The catalogs 
cover a cubic volume with side 1.5/i _1 Gpc. To construct 
an "observed" catalog, we put the observer in one corner 
of the box and consider the subsample of galaxies in the 
shell octant between the observer's z = 0.5 — 0.55. For 
simplicity, we do not apply photometric redshift errors 
nor do we introduce redshift shifts due to peculiar veloc- 
ities. This latter effect would only be significant on very 
large scales anyway (see Figure [2]). Each mock covers 
7r/2 rad 2 and consists of about 125,000 galaxies. Since 
both area and galaxy number are thus roughly half the 
values for the z = 0.5 — 0.55 redshift bin of the true 
data, the number density is comparable to that of the 
true photometric sample. We apply this procedure for 
eight different corners to get eight different lines of sight 
per simulation. Note, however, that these lines of sight 
are not completely independent as they are based on the 
same simulation volume. Finally, for each line of sight, 
we average the spectra over all twenty independent real- 
izations in order to increase the signal to noise ratio. The 
covariance matrix for the mock angular power spectrum 
is rescaled accordingly to reflect the decrease in covari- 
ance due to taking the average. 

As described in section [3j we consider a model char- 
acterized by the cosmological parameters and a galaxy 
bias bo (giving our mock bin the label 0), and a more 
conservative model with bias b and nuisance parameter 
do, so that the spectrum is given by 



C, 



(00) 



k 2 dkP m {k,z = 0) (Aj; l) (fc)) 2 + a 



(15) 

In the galaxy bias-only version, clq is simply set to zero. 

As a direct test of this model, we fit it to the aver- 
aged mock spectrum. In this first approach, we keep the 
cosmology fixed to the mock's input cosmology and re- 
strict the fit to the range I — 30 — 200. The only free 
parameters are thus either bo or (&o,ao)- We use a mod- 
ification of the publicly available COSMOMC package 



Lewis & Bridle (2002) to sample this parameter space 
using Monte Carlo Markov Chains (MCMC). We show 
the resulting best-fit spectra together with the mock av- 
erage in Fig. [5] Considering first the default, galaxy bias- 
only model (Black curve) , we find that the best fit to the 
mock result has a linear bias 60 = 2.02 (with uncertainty 
a(b ) < 0.01) and has \ 2 = H-3. This should be com- 
pared to an expectation value of (x 2 ) — 16 based on 17 
data bins and one free parameter. The galaxy bias model 
thus provides a good fit to the simulated spectrum (the 
probability of getting a \ 2 lower than 11.3 for an expec- 
tation value of 16 is approximately 20%). 

Next, including the shot noise- like parameter ao to 
take into account potential residual shot noise and/or 
non-linear effects not captured by our simple Halofit plus 
scale-independent galaxy bias model, we find a best fit 
model with ao = 1.1 • 10~ 6 and 60 = 2.00. However, the 
uncertainty in ao is c(ao) = 1.0 • 10~ 6 so the preference 
for a non-zero value cannot be considered significant. In 
this model, we find \ 2 = 10.1, to be compared to the 
expected (\ 2 ) = 15. This is only a marginal (Ax 2 ~ 1.2) 
improvement. 

Restricting the fitting range to I — 30 — 150, we find 
Ax 2 — 1-1 between the two best-fit models, and ao = 
(1.8 ± 1.9) x 10~ 6 . 

The comparison above suggests that for the range 
£ = 30 — 200, the galaxy bias model without an extra 
nuisance parameter may be sufficient. We now undertake 
a more complete check of our model and the entire cos- 
mology analysis by using MCMC to fit the full space of 
cosmological parameters as well as the galaxy bias (and 
shot noise parameter) to the averaged mock spectrum. 
The differences between the resulting best-fit parameter 
values and the "true", i.e. input, values give an indica- 
tion of the parameter bias introduced by our method. 
To break parameter degeneracies, while not letting the 
prior bias us away from the input cosmology, we include 
a "mock" CMB prioip^| which will provide a likelihood 
similar to the true WMAP7 one, except shifted to be 
centered around the mock input parameters. 

We want any deviation between the input cosmology 
and the recovered cosmology to be small compared to the 
uncertainty level of the actual data. We therefore take 
information from our final results and use the uncertain- 
ties for the WMAP+HST+CMASS (4 iax = 200) real 
data case for comparison. If the biases on parameter esti- 
mation are small compared to these numbers, it provides 
strong motivation for considering our approach sound, as 
most uncertainties in the next section will be larger than 
in the WMAP+HST+CMASS case. Therefore, param- 
eter uncertainties a referred to in the remainder of this 
section are these data-based uncertainties. 

Starting with the parameter spacj 2 ^] 
{flbh 2 , fioM/i 2 , 0, A s , n s , t, 6 }, we find that all cos- 
mological parameters are reproduced to within la of 

19 The "mock" CMB prior is defined by x|vMAP7 — (P« — 
Pi, input) COV" 1 (pj — p.input)) where pi are the parameters at each 
point of the chain, p; ; i npu t the input parameters, COVij is the 
WMAP7 covariance matrix and i,j are summed over. 

20 The parameters 8, A s and r represent the ratio between the 
sound horizon and the angular diameter distance at decoupling, 
the scalar amplitude of primordial fluctuations and the reionization 
optical depth, respectively. 
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Fig. 4. — Effect of neutrinos on the angular power spectra. The solid and dashed curves depict the massless and Em„ 
respectively. 



0.3 eV cases, 



the input values (although the parameter most affected 
by the mock CMASS data, fl^Mh 2 , is higher than the 
input by close to la). 

Unfortunately, we do not have mocks based on 
a cosmology with non-zero £m„. One check we 
can do, however, is to fit a model with parameters 
{ilbh 2 ,H,DMh 2 ,9, A s ,n s ,T,T,m u ,bo} to our T,m v = 
mock spectra. The parameters affected by far the most 
by the angular spectra are (again) J7dm^ 2 and Em„. We 
show the posteriors of this calculation in Fig. [61 In the 



left panel, the vertical lines indicate the Odm^ input 
value, and the la and 2a bounds based on the uncer- 
tainty a from the actual data. The blue points with 
error bars are the posterior mean values and la recov- 
ered errors after fitting to the averaged mock spectrum. 
Note that the recovered error bars (from the averaged 
mock power spectrum) are typically similar to the data- 
based error bars. While the different lines of sight are 
not entirely independent, Fig. [6] points towards a bias of 
about 1 — 1.5a in firjM^ 2 - F° r the neutrino mass, the 
right panel shows the posterior probability distributions 
in blue. The posteriors are always consistent with the in- 
put value £m„ = and can be interpreted as providing 
upper bounds. We have made the same plot as in the 
left panel for the other parameters and they were biased 



significantly less (as their reconstruction is dominated by 
the mock CMB prior). 

Adding the nuisance parameter an, we obtain the red 
points and curves in Fig. [6] The effect of marginalizing 
over an is to diminish the parameter bias so that f^DM^ 1 
is typically reconstructed to well within la. We at- 
tribute this change to an accounting for a possible scale- 
dependence in galaxy bias on quasilinear scales. The 
neutrino constraints are also still consistent with the in- 
put, although the mock upper limits do become signif- 
icantly weaker. We have also studied mock cosmology 
constraints using £ max = 150, and found that the main ef- 
fect is to widen the posterior distributions slightly, while 
the change in parameter bias relative to f max = 200 is 
small. 

We conclude that our galaxy bias-only model and the 
fitting method used here properly reproduce the input 
cosmology for our choices of £ max , except that there is 
a bias of about 1 — 1.5a in f^DM^ 2 - The model with 
nuisance parameter an removes this parameter bias at 
the cost of larger error bars. While the bias in 51dm^ 2 
is not extreme, being only slightly above the la level, it 
is sufficiently worrying that we will quote results for the 
galaxy bias-only model and for the more conservative 
model with shot noise-like parameters. Changing £ ma , x 
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Fig. 5. — Example of an averaged mock spectrum (green points with error bars) and theoretical spectra (solid lines). Fixing the cosmology 
to the mock input cosmology (see text), we fit the averaged mock spectrum in the range I = 30 — 200 to our model described in the text. 
The black curve is the resulting best-fit spectrum if we only fit a (scale-independent) galaxy bias bo (best-fit value bo = 2.02), while the 
red curve is the best fit in a model that also includes the nuisance parameter ao (best fit values &o = 2.00, ao = 1 .05 ■ f0 — 6 ). To provide a 
hint of the importance of non-linear effects in this multipole range, we plot the spectrum based on a linear three dimensional matter power 
spectrum in blue (bo = 2.02, ao = 0) 

between 150 and 200 does not have a large effect on how 
well the models compare to mocks, suggesting that both 
are reasonable choices. We will quote results for both 
ranges. 

6. RESULTS 



pilation Amanullah et al ( |2010 ) and the meas urement of 
the BAO scale based on SDSS Data Release 7|Abazaiian 



et al. ( 2009 ) spectroscopic data from Ref . Percival et al 



While the CMASS angular galaxy power spectra carry 
useful information about the sum of neutrino masses 
Em„, the effect of Y>m v is degenerate with certain other 
parameters which are not well constrained from the an- 
gular spectra alone. There are many combinations of 
external data sets that our angular spectra can be com- 
bined with to fix this problem. One approach would be to 
optimize the neutrino bound by combining as many data 
sets as possible. However, we choose instead to focus 
as much as possible on the effect of the CMASS photo- 
metric data and therefore consider mostly simple p riors. 
Our two main pr ior choices are WMAP7 CMB data |Lar^] 
son et al.| ( |2010[ ) and the combination of WMAP 7 with 
the MSI m easurement of the Hubble parameter |Riess| 
et al. (20111. At the end of this section, we will briefly 
consider the effect of adding the Union 2 supernova com- 



We again use a modif ication of the public ly avail- 
able COSMOMC package [Lewis fc Bridle| p002| ) to sam- 
ple the parameter space using MCJVIO TJur parame- 
ter space consists of the six usual ACDM parameters, 
(i7f,/i 2 , f^DM^ 2 , 0, ln(10 10 A s ), n s , r), the neutrino mass 
fraction /„, defined as fl^/fl^M (where SJdm includes 
cold dark matter and massive neutrinos), in addition 
to Asz, describing the amplitude relative to a tem plate 
of the Sunyaev-Zel'doy ich contribution to the CMB |Ko- 



matsu & Seljak (2002), the four galaxy bias parameters 
bi and (optionally) the four nuisance parameters , leav- 
ing us with a maximum total number of parameters of 
sixteen parameters. We put uniform priors on these pa- 
rameters and derive Em„ using Eq. (IT]). 

We first consider the WMAP7 prior and show how 
the neutrino bound improves as CMASS data are added. 
The resulting 95% CL upper limits are shown in the top 



10 




Fig. 6. — Left panel: Recovered values of Odm^ 2 from averaged mock spectrum together with CMB prior. We consider spectra from eight 
different lines of sight. The points with error bars show the posterior mean values and la error bars after the Monte Carlo analysis for two 
scenarios: varying £m„ with (red) and without (blue) ao marginalized. The vertical magenta lines indicate the input value SIfjm' 1 = 0.11166 
(solid) and the input ± one and two a, where a = 0.0032 is the parameter uncertainty based on the data set WMAP7+HST+angular 
spectra (a^ = fixed). There is a bias of about la without the nuisance parameter, which disappears when a; is marginalized over. Right 
panel: The posterior neutrino mass distributions for the two cases discussed above. The mock constraints are consistent with the input 
cosmology of Em„ = 0. Other parameters are all reconstructed to close to their input values and are not strongly affected by the angular 
spectra. 



95% CL £m„[eV] 


prior only 


prior+CMASS,4mx = 150 


prior+CMASS/ max = 200 


WMAP7 prior 


1.1 


0.74 (0.92) 


0.56 (0.90) 


WMAP7 + HST prior 


0.44 


0.31 (0.40) 


0.26 (0.36) 


TABLE 1 



The 95% confidence level upper limits on the sum of the neutrino masses Sm„. The top row investigates the effect of 

ADDING THE CMASS GALAXY POWER SPECTRA TO A WMAP PRIOR WHILE THE BOTTOM ROW USES WMAP AND THE Hq CONSTRAINT 
FROM HST AS A PRIOR. IN PARENTHESES WE SHOW RESULTS FOR THE MORE CONSERVATIVE MODEL MARGINALIZING OVER THE SHOT 

NOISE-LIKE PARAMETERS 0^. 



row of Table [TJ with the results with <Zj marginalized in 
parentheses. The bound improves from 1.1 eV for CMB 
only to 0.56 eV for CMB with CMASS data (£ max = 200). 
This constrain t is comparable t o the limit £m„ < 0.62 
eV derived by Reid et al. ( 2010[ ) from the DR7 spectro- 
scopic sample/ It thus appears that the advantage of 
spectroscopic redshifts (providing information on clus- 
tering along the line of sight) in that sample is offset 
by the advantage of the current sample having a larger 
volume, although there are other differences between the 
samples and analyses as well. Note, however, that the 
constraint deteriorates significantly when marginalizing 
over the nuisance parameters eij. In this case, the mass 
bound is not significantly better than with CMB alone. 
We show the posterior probability distributions for Em„ 
and the other cosmological parameters in Fig. [7[ 

We next consider the constraints using WMAP7 with 
HST Hq prior. The CMB alone provides a strong mea- 
surement of one combination of late-universe parameters 
through its sensitivity to the distance to the last scatter- 
ing surface. However, this distance measurement leaves 
a degeneracy between fl\ and Sm^ so that the CMB- 
only limit on the neutrino mass arises mainly from the 
effect of neutrinos on the primary anisotropics and not 
from this distance measurement. Measuring H con- 
strains a different combination of late universe param- 
eters and thus breaks the CMB degeneracy. This is why 



the WMAP7+HST bound is so much stronger than the 
WMAP7-only one, i.e. £to„ < 0.44 eV as opposed to 
Y>m v < 1.1 eV. Adding the CMASS angular spectra 
tightens the bound significantly so that an impressive up- 
per bound of Y>m u < 0.26 eV is reached for £ max = 200 
(in the bias-only model), as is shown in the second row 
of Table [Tj The effect of marginalization over dj is again 
to bring the constraint back to closer to the CMB+HST 
bound. 

The posteriors for all cosmological parameters are 
shown in Fig. [SJ In addition to the full likelihoods 
for Em„, summarized in Table [TJ the JJdm'i 2 posteri- 
ors are worth noting. The effect of the angular spec- 
tra is to strongly shift the average value of this parame- 
ter (blue curve) , while including the nuisance parameters 
(red curve) weakens the shift. The last two parameters in 
Fig. [8] (and [7|, Oa and H 0l are not independent and can 
be expressea in terms of the preceding parameters. The 
shift thus is really only significant for one independent 
parameter, J7dm^- 2 , in our basis. The results in section 
[3] suggest that the shift in the bias-only case might par- 
tially be a bias due to our model and that the results 
with a,i marginalized are unbiased. 

We do not explicitly show the correlations between pa- 
rameters, but have verified that, in the CMB+CMASS 
case, the neutrino mass has its strongest degeneracies 
with fiDM^ 2 , the bias parameters bi and erg. While, in 
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Fig. 7. — Cosmological constraints with a WMAP7 CMB prior. We show the probability distribution functions for CMB only (black), 
CMB with CMASS spectra in the range I = 30 - 150 (blue dashed) and CMB with CMASS spectra in the range t = 30 - 200 (blue solid). 
The red curves represent the constraints in the conservative model where we marginalize over a set of nuisance parameters a^. 

agreement with our discussion in sectionHJ the inclusion 
of the Hubble prior removes the £m„ — SI^m^ 2 degener- 
acy, the strong correlations with hi and erg remain. 

We have also added supernova and BAO data to the 
CMB+HST+CMASS data set, and considered the neu- 
trino mass bound in the bias-only model, but we found 
negligible improvement (from 0.26 eV to 0.25 eV) rela- 
tive to the case without these additional data sets. These 
additional data sets do carry significant information, but 
this information is degenerate with the information al- 
ready present in the three default data sets. 

For the multipole range £ = 30 — 150, we show the 
results using dashed lines in Figures [7] and [8) The 95% 
CL upper limit for CMB+HST+CMASS in this case is 
0.31(0.40) eV and for CMB+CMASS it is 0.74(0.92) eV 
fixing (varying) a, = 0. A significant amount of informa- 
tion is thus contained in the large I range of multipoles, 
which makes sense as the number of modes is large. 

Finally, we consider the question of where most of 
the neutrino mass information comes from. In principle, 
massive neutrinos affect the angular power spectra both 
by their small-scale suppression of the three-dimensional 
power spectrum, and by changing the projection of phys- 
ical scales onto angular scales through their effect on the 



background expansion. As discussed in section |4j we ex- 
pect the former effect to carry more constraining power 
than the latter effect. We test this by running Monte 
Carlo chains where the effect of massive neutrinos on the 
three-dimensional power spectrum is artificially taken 
out, while the effect on the background expansion is left 
intact. Specifically, we replace the usual linear CAMB 
power spectrum by the spectru m given by the Eisenstein 
and Hu (EH) fitting formula lEisenstein fc Hu| ( |1998[ ), 
which does not include the effect of massive neutrinos. 
We find that in this setup, including the CMASS galaxy 
power spectra does not improve the neutrino mass bound 
relative to the case with CMB, or with CMB+HST, only. 
In other words, the projection effect alone carries little 
information on neutrino mass (at least after marginaliz- 
ing over the effects of other parameters) and the bounds 
quoted in this manuscript can be attributed to the small- 
scale suppression of the three-dimensional power spec- 
trum. 

7. CONCLUSIONS 

We have exploited angular power spectra from the 
SDSS-III DR8 sample of photometric galaxies with 
CMASS selection criteria to put interesting constraints 
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Fig. 8. — Cosmological constraints with a WMAP7 CMB and Hubble parameter prior. We show the probability distribution functions 
for CMB + H only (black), CMB + H with CMASS spectra in the range £ = 30 - 150 (blue dashed) and CMB + Ho with CMASS 
spectra in the range i = 30 — 200 (blue solid). The red curves represent the constraints in the conservative model where we marginalize 
over a set of nuisance parameters a^. While this marginalization degrades the neutrino bound, simulations have shown it removes the bias 
in UoMh 2 (see section [B|. 



on the sum of neutrino masses. We have used mock 
galaxy catalogs based on N-body simulations and HOD 
modeling to test two models for the angular galaxy spec- 
tra. Based on these tests, we decided to compare the 
data to theoretical spectra based on the non-linear mat- 
ter power spectrum augmented by a linear galaxy bias 
factor. However, since this model does result in a bias 
in Adm/» 2 of ~ 1 — 1.5a, we have also fitted the data 
to a more conservative model, with an additional set of 
shot noise-like fitting parameters, in which this bias is 
virtually absent. The tests also motivated us to use the 
multipole range I = 30 — 200, but we quoted results for 
the more conservative choice I = 30 — 150 as well. The 
added advantage is that this analysis provides insight 
into the range of scales that yields the galaxy clustering 
information. 

Combining the CMASS data with a CMB prior from 
the WMAP7 survey, we find an upper bound £m„ < 0.56 
eV (0.90 eV) at 95% confidence level for £ max = 200 in 
the model with free parameters hi {bi and a{). Adding the 
HST measurement of the Hubble parameter, the proba- 
bility distribution tightens and we find £m„ < 0.26 eV 



(0.36 eV). We have also considered the effect of adding 
supernova and a (lower redshift) BAO measurement, but 
when the HST prior is included already, these addi- 
tions lower the upper limit to 0.25 eV (in the bias-only 
model). Considering the dependence on the multipole 
range, characterized by a maximum multipole i max , we 
find that a significant amount of information resides in 
the largest multipoles £ = 150 — 200, but that even for 
4m = 150, the galaxy spectra place a strong bound 
on neutrino mass. Our main results are summarized in 
Table Q] 

It is interesting to compare these results to the outcome 
of an analysis of a similar (but smaller) high redshift 
SPSS photom etric catalog, the MegaZ s ample Collister 



et al.| (|2007j). In |Thomas et al.| (|2010|, the strongest 
bound quoted is a 957c CL upper limit of 0.28 eV, in- 
cluding SN and BAO data in addition to CMB, HST 
and MegaZ. However, this particular bound is based on 
a multipole range with £ max = 300 and no nuisance pa- 
rameters dj . As we have discussed extensively, we believe 
4m = 200 (or even slightly lower) is a better choice if 
one wants to avoid significant, unknown non- linear cor- 
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rections to the galaxy bias. For this choice, the MegaZ 
sample yields an upper bound of 0.34 eV. Assuming the 
aggressive, bias-only model (as in the MegaZ analysis), 
the value we find for the CMASS sample is 0.25 eV, which 
is thus a significant improvement. However, it must be 
kept in mind that this model causes a small parame- 
ter bias and that the more conservative model yields a 
weaker bound of 0.36 eV. 

The bounds presented here rule out the quasi- 
degenerate neutrino mass hierarchy. For example, for 
Zm u = 0.25 eV, it follows from |Am| 3 | = 2.5 • 10~ 3 eV 2 
that the mass difference |m 3 — ma| ~ 0.015 eV, so that 
the largest mass difference is \rri3 — m2|/(£m„/3) « 20% 
of the average neutrino mass. We are thus entering the 
regime where the mass splittings are significant. Looking 
forward, the prospects are exciting. As the sensitivity of 
cosmological neutrino mass measurements improves, the 
sum of the masses will either be measured, i.e. a value 
that can be distinguished from zero will be found, or the 
upper limit will be sharpened. However, even in the lat- 
ter case, interesting things can be learned. If the sum of 
the masses is found to be less than ~ O.leV, this rules 
out the inverted hierarchy, leaving the normal hierarchy 
as the only option. Moreover, even in the normal hier- 
archy, Em„ is not allowed to be lower than ~ 0.05eV 
so that sooner or later a measurement, rather than an 
upper bound, can be expected. 

Finally, we note that BOSS is currently taking spectra 
for a sample of galaxies with the same selection criteria 
as the galaxies considered in this paper and will reach 
a sample of approximately 1, 500, 000 galaxies when fin- 
ished in 2014. The results presented here and in our 
companion papers are thus only the tip of the iceberg of 
what can be done with BOSS. The spectroscopic data 
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